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We investigate the QCD phase diagram in the strong coupling limit by using a newly developed 
auxiliary field Monte-Carlo (AFMC) method. Starting from an effective action in the leading 
order of the 1 /g 2 and 1 jd expansion with one species of unrooted staggered fermion, we solve the 
many-body problem exactly by introducing the auxiliary fields and integrating out the temporal 
links and quark fields. We have a sign problem in AFMC, which is different from the original 
one in finite density lattice QCD. For low momentum auxiliary field modes, a complex phase 
cancellation mechanism exists, and the sign problem is not serious on a small lattice. Compared 
with the mean field results, the transition temperature is found to be reduced by around 10 % and 
the hadron phase is found to be extended in the larger chemical potential direction by around 20 
%, as observed in the monomer-dimer-polymer (MDP) simulations. 
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1. Introduction 

QCD phase diagram is expected to have rich structure, extensively studied in current high- 
energy heavy ion collision experiments, and closely related to compact astrophysical objects and 
phenomena |jl|]. In the early universe as well as in heavy ion collisions at collider energies, the 
quark gluon plasma evolves to hadronic matter via the crossover transition at small quark chemical 
potential /I and high temperature T. The first order phase transition from hadronic to quark matter 
may be realized in cold dense matter such as the neutron star core. The QCD critical point (CP) 
connects these crossover and first order transitions, and characteristic large fluctuations of the order 
parameter may be observed in heavy-ion collisions or in black hole formation processes [gj. 

The phase transition at finite density is, unfortunately, much less known compared with the 
low density crossover transition. The sign problem of lattice QCD at finite jU makes it difficult to 
perform precise calculations of dense cold matter. The cancellation of the statistical weight given 
by the fermion determinant may not be just a technical problem, as recent studies suggest that we 
cannot reach CP in phase quenched simulations [Qj. Therefore we need to find other methods than 
the phase quenched simulation in order to directly sample appropriate configurations in dense cold 
matter for the discussion of CP and the first order transition. 

One of the hopes may be found in the strong coupling lattice QCD. In the strong coupling 
region, we can carry out the link variable integral prior to the fermion integral for each order of 
the inverse squared coupling l/g 2 The QCD phase diagram has been studied under the mean 
field approximation in the strong coupling limit (leading order, 0(l/g )) [||], the next-to-leading 
order (NLO, 0(1/ g 2 )) [§], and the next-to-next-to-leading order (NNLO, 0(l/g 4 )) [0]. We can 
also take account of the fluctuation effects beyond the mean field approximation. We obtain the 
effective action of quarks after the link variable integral, and the fermion integral leads to the 
partition function in the form of sum over monomer-dimer-polymer (MDP) configurations [j|]. 
The fluctuation effects are found to modify the phase diagram shape moderately: The transition 
temperature is shifted to lower T direction by 10-20 %, and the hadronic phase expands to higher 
/i direction by 20-30 % M\. Until now, MDP simulation has been performed only in the strong 
coupling limit, 1 /g 2 = 0. Since both finite coupling and fluctuation effects are important to discuss 
the QCD phase diagram, we need to develop a theoretical framework which includes both of these 
effects. 

In this work, we develop an auxiliary field Monte-Carlo (AFMC) method for the strong cou- 
pling lattice QCD. In AFMC, we carry out the integral over the auxiliary fields, which are intro- 
duced to decompose the fermion composite product. In the mean field method, auxiliary fields are 
introduced and assumed to be static and constant. Thus AFMC is a straightforward extension of 
the mean field method, and may be also applicable to finite coupling cases. We discuss here the 
QCD phase diagram in the strong coupling limit. 

2. Auxiliary field Monte-Carlo method 

We consider an asymmetric lattice {a % =a/y) of size L 3 x N x for color SU(Af c = 3) with one 
species of unrooted staggered fermion in the strong coupling limit (SCL). Throughout this paper, 
we work in the lattice unit a = 1 , where a is the spatial lattice spacing. The lattice QCD action is 
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given as 



>LQCD 



V: 



7 



(2.1) 



(2.2) 



where £ r and U VjX represent the quark field and the link variable, respectively, rij x = (— l) VoH KVM 
is the staggered sign factor, and chemical potential pL is introduced in the form of the temporal 
component of a vector potential. In SCL, we can ignore the plaquette action terms, which are pro- 
portional to l/g 2 . By integrating out spatial link variables, we obtain the SCL effective action [|H]]. 



>eff 



(2.3) 



Here we adopt the effective action in the leading order of the 1 jd expansion, where d is the spatial 
dimension, d = 3. The nearest neighbor four fermi interaction, the second term in Eq. (2.3), is 
rewritten in the momentum representation as, 

1 ^ L 3 



4N c r 2 % 



4N c f 



£/(k)M_ k , t M M , /(k) = £cos kj 



(2.4) 



k.r 



where the Fourier transformation is defined as M x =r x )T ) = Lk e ' kx ^k,T- 

We shall now introduce auxiliary fields and decompose the interaction terms. We apply the ex- 
tended Hubbard-Stratonovich (EHS) (or Hubbard-Stratonovich-Miura-Nakano-Ohnishi-Kawamoto) 
transformation [Q], 

e aAB = J dwdw * e -a\rW-AW-rB\ ; q.5) 

where dxj/dxj/* = JRe^JImi//. By introducing two auxiliary fields simultaneously, we canbosonize 
any kind of composite product. For the interaction term Eq. (2.4), we find that the positive and neg- 
ative meson hopping matrix eigenvalues appear in pair, /(£) = — /(k) with k = k + (jt, 71, n). We 
refer to the modes with momentum k satisfying /(k) > and /(k) < as positive and negative 
modes, respectively. The bosonization of these terms is carried out as 

exp (o/(k) [M_ kiT M kiT -M_ £iT M £jT ] } 

= J dak^da^d^d^expl-afik) [\<y k .t\ 2 + \K k . t \ 2 

+a k * T M kiT +M_ k , T a k , T -/(-) T 7r^ T M kT -/(-) T M_ kiT 7r k , T ]} , (2.6) 

where k is chosen to be one of the positive modes, /(k) > 0. By construction, a kT and 7T k T satisfy 
the relation a kT = c k T and 7T_ k T = T . The bosonized interaction terms are given as 

ST = \ I W - V-] +£m x M x + 



4N c y 2 



I /(k) 

k,T,/(k)>0 



|<7 k . T | 2 + |^k.r| 2 



ffl. 



1 



7 4N c y 2 



£ (a + /£7r) x+; + (a + /£7r). v _; 



°x= 

k,/(k)>0 



k,/(k)>0 



(2.7) 

(2.8) 
(2.9) 
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where e x = (— )*o+*i+*2+*3 corresponds to ^5 = 75 <g> 75 in the spinor-taste space. The lattice QCD 
action Eq. ( jLl| ) is invariant under the chiral U(l) transformation, % x — > e ie * G % x , and the chiral 
transformation mixes Ok T and n^.z- Thus Ok T and % T at small k are regarded as the usual chiral 
(a) and Nambu-Goldstone (71) fields, respectively. 

We can carry out the Grassmann and temporal link (Uq) integrals semi-analytically. 



r-AF 
°eff 



M,/(k)>0 



L 3 m 

AN c y 2 



£log [X Nt (x) 3 -2X Ni (x) + 2cosh(3N^/r 2 )] , 

X 

(2.10) 



where X^ t {x) is a known function of m x [11], and we can obtain it by using a recursion formula 
with N T steps. When m x -t x ^) is independent of T, we obtain X^ T = 2cosh(A r T arcsinh m*). 

We perform Monte-Carlo integral calculations over the auxiliary fields (Ck.f, %,t) based on 



the auxiliary field effective action Eq. ( |2.10| ). We refer to this treatment as an auxiliary field Monte - 



Carlo (AFMC) method. We have made two approximations to obtain the effective action Eq. ( |2.3| ) 
(leading orders of strong coupling expansion and l/d expansion), whereas no approximations have 



been invoked to calculate observables based on 5 e ff in Eq. fl2.3|). One of the merit of AFMC in the 
strong coupling limit is that the fermion matrix is decomposed into that at each spatial site. The 
numerical cost is proportional to the space-time lattice volume multiplied by the one-dimensional 
size, L 3 xN t x (L or N t ). As in fermion many-body problems in other fields of physics, unfortu- 
nately, we have a sign problem: The fermion self-energy m x is complex, and the second term in the 



auxiliary field effective action Eq. (2. 10) contains the imaginary part. As a result, the statistical 



weight exp(— ) has a phase, coming from the negative modes % This sign problem is weak- 
ened in part by the phase cancellation mechanism. Since negative modes involve ie x , the phase on 
one site from low momentum 7i k T modes is tend to be cancelled by the phase on the nearest neigh- 
bor site. Phase from high momentum modes does not cancel, but we expect that high momentum 
modes are less relevant to long wave phenomena such as the phase transition. Nevertheless, we 
demonstrate that AFMC works in a small lattice such as 4 3 x and 6 3 x N T . 



3. Phase diagram in AFMC 

We discuss here the AFMC results of the phase diagram in the chiral limit (wiq = 0) on a small 



lattice 4 3 x N x and 6 3 xiV T . Following the arguments in [12], we assume the temperature is given 
as T = fjN x . 

In the left panel of Fig. [j], we show the average sign factor (cos0), the root mean square 
chiral condensate = ^/( a k=o z + K k=o t)> an( ^ tne 1 uar k number density p q , as functions of the 
temperature on a 4 3 x 4 lattice. Lines connect fixed fi/T results. At small \ijT, (j) and p q smoothly 
change around the transition temperature T c . During this change, the ty 1 distribution is always 
single peaked as shown in the right upper panel of Fig. [T[ The peak at a finite (j) at lower T moves 
to the peak at = at T c . This behavior suggests that this transition is a would-be second order 
transition, which could be the second order for a large enough volume. At large jl/T, we find 
differences between the results starting from the Wigner phase and the Nambu-Goldstone phase 
initial conditions (dashed lines), as far as the sampling steps are not very large. In the initial 
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conditions, we set the auxiliary fields static and constant. We set o x = 0.01 and 1 in the Wigner 
and Nambu-Goldstone phase initial conditions, respectively, and % x = is taken in both types of 
initial conditions. The 2 distribution becomes double peaked as shown in the right bottom panel 
of Fig. p]. The initial condition dependence suggests the two states are separated by a high barrier 
in the effective potential, and the transition is suggested to be would-be first order. 

The average sign factor is large enough, (cos 6) > 0.9, and there is practically no sign problem 
on a small lattice. Generally, the average sign factor increases with increasing T for a fixed fugac- 
ity, and increases with increasing jl except for the transition region. The global /I dependence is 
understood from the form of the effective action. The imaginary part in comes from the auxil- 
iary fields, and their effects becomes smaller at larger jj, as we can find in Eq. ( 2.1C ). Around the 
transition region, finite momentum modes will contribute, where the phase cancellation mechanism 
does not work completely, and the average sign factor is suppressed. 

We have evaluated the transition temperature based on the chiral susceptibility peak in the 
second order region of fl/T, and by comparing the effective potential value at the local minima in 
the first order region. In Fig. |2[ we compare the AFMC phase boundary with that in the mean field 
approximation [J|] and in the MDP simulation As discussed in the previous works using the 
MDP simulation the transition temperature is found to be reduced at jU ~ by around 10 %, 
and the hadron phase is found to be extended in the larger /I direction at low T. 

In the strong coupling limit, the (would-be) first order phase boundary is insensitive to the 
spatial lattice size; the phase boundary on the 6 3 x 4 lattice is almost the same as that on the 4 3 x 4 




Figure 1: Average sign factor (top), chiral condensate (middle), and quark number density (bottom) as 
functions of temperature on a 4 4 lattice. 
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SCL (l/g 2 =0) 
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Figure 2: Phase diagram 

lattice, while the transition is found to be shaper on the 6 3 x 4 lattice. For larger N T lattices, the 
transition chemical potential is found to be larger. Extrapolated results to N r = °° is consistent 
with the continuous time MDP simulation [0]. These observations imply that AFMC can be a 
promising tool to discuss finite density lattice QCD, and it can be an alternative method of the 
MDP simulation. 

4. Summary 

We have discussed the QCD phase diagram by using the auxiliary field Monte-Carlo (AFMC) 
method. Starting from an effective action 5 e ff (X>Xi ^o) of quark and temporal link variables in the 
leading order of the l/g 2 and 1/d expansion with one species of unrooted staggered fermion, we 
obtain an effective action S™ s (x,X, Uo, Ck jT , 7Tk.-r) in the bi-linear form of quarks via the extended 
Hubbard-Stratonovich transformation. Integral over the quark and temporal link variables is car- 
ried out analytically, and the auxiliary field effective action 5^(ak,T 5 %.t) is obtained. Finally, we 
integrate out the auxiliary fields (Ok,T 5 %.t) using the Monte-Carlo method. This procedure corre- 
sponds to solving the many-body problem of 5 e g exactly. While we have a sign problem in AFMC, 
the problem is not serious on a small lattice partly due to the phase cancellation mechanism of the 
low momentum % T modes. The phase diagram in AFMC is found to be compatible with that in 
the monomer-dimer-polymer (MDP) simulations [^]. Thus the QCD phase diagram in the strong 
coupling limit in MDP is confirmed to be correct. 

We can extend AFMC to include finite coupling effects in a straightforward manner. The sign 
problem with finite coupling or on a larger lattice, however, may be more severe. The present sign 
problem comes from high momentum negative modes, then it is an interesting direction to study to 
integrate out high momentum auxiliary field modes above a given cutoff in an approximate manner 
such as in the saddle point method. Since the auxiliary fields introduced here are color singlet, the 
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momentum cutoff does not violate the gauge symmetry. If the cutoff dependence is small, we hope 
that it would be possible to evaluate the phase diagram on a large lattice at finite coupling with less 
severe sign problem. 
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